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Abstract 

Quantum field theory on non-commutative spaces does not enjoy the usual ultraviolet- 
infrared decoupling that forms the basis for conventional renormalization. The high 
momentum contributions to loop integrations can lead to unfamiliar long distance 
behavior which can potentially undermine naive expectations for the IR behavior of 
the theory. These "anomalies" involve non-analytic behavior in the noncommuta- 
tivity parameter 9 making the limit 6^0 singular. 

In this paper we will analyze such effects in the one loop approximation to gauge 
theories on non-commutative space. We will see that contrary to expectations poles 
in 6 do occur and lead to large discrepancies between the expected and actual infrared 
behavior. We find that poles in 9 are absent in supersymmetric theories. The 
"anomalies" are generally still present, but only at the logarithmic level. A notable 
exception is non-commutative super Yang Mills theory with 16 real supercharges in 
which anomalous effects seem to be absent altogether. 
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1 The UV/IR Connection in Non-commutative Field 
Theory 



Field theories formulated on non-commutative spaces are interesting in both their own 
right as well as for their applications to string and matrix theories PJH. These theories 
are characterized by a non-commutativity parameter 9 with dimensions of length squared. 
Classically and in the tree-level approximation the behavior of the theory for momenta 
much less than Q~ x l 2 is the same as for the corresponding commutative theory. However 
this is not necessarily the case for the quantum theory. The non-commutativity can lead to 
unfamiliar effects of the ultraviolet modes on the infrared behavior which have no analog 
in conventional quantum field theory ||. 

The origin of the strange mixing of IR and UV effects in non-commutative field theory 
can be understood in a simple way [|J Q. The field quanta in such a theory can be thought 
of as pairs of opposite charges moving in a strong magnetic field . The spatial locations 
of the two charges are defined by a center of mass position x l cm and a relative coordinate 
A m . The relative coordinate is related to the spatial momentum p by 



where 6 is an antisymmetric matrix with components in the spatial directions. In this 
paper we will consider the case of 3 dimensional space. Without loss of generality 9 can 
be taken to lie in the (1,2) plane 



The momentum in the (1, 2) plane will be called P. We will also use the notation Pi = 
OijPi . Thus a particle moving with momentum P along the X 1 axis has a spatial extension 
of size \9P\ in the X 2 direction. The growth of the size of a particle with its momentum 
has interesting consequences. For example, when a quantum of momentum P scatters off 
a target at rest, the scattering amplitude will spread in impact parameter space over a 
distance \6P\. 

There are also important and somewhat bizarre consequences for Feynman loop inte- 
grations. Roughly speaking, when a particle of momentum P circulates in a loop it can 
induce an effect at distance \6P\. The high momentum end of the integrals can give rise to 
power law long range forces which are absent entirely in the classical theory. We will call 



(i.i) 




(1.2) 
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Figure 1: Mass renormalization correction in 4 thoery. 

such effects "anomalies" although we emphasize that they do not signal an inconsistency 
in the theory but rather a violation of naive expectations. 

In the case of planar diagrams the effect of non-commutativity is simple. Every diagram 
gets multiplied by a phase factor that depends only on the momenta of the external lines. 
Thus the Feynman integrals are exactly as in the commutative theory. In the non-planar 
case the situation is more interesting. The phase factors now involve products of external 
momenta p and internal momenta I in the form 

e ipm = e ipI . (1.3) 

If the diagram in question is divergent in the commutative theory, the effect of the oscil- 
lating phases is typically to regulate the diagram and render it finite. But as P — > the 
phases become ineffective and the diagram diverges at P = 0. This is the mechanism de- 
scribed in detail in ||. We will begin by reviewing a simple example from non- commutative 
</> 4 theory. 

The diagram in question, fig(l), is the lowest order mass renormalization correction to 
the propagator. We are interested in the non-planar contribution which in the commutative 
theory has the form 

The diagram is quadratically divergent and is renormalized by a mass counter term. 

In the non-commutative case the integrand has an additional factor exp(ipl) where p 
and I are the external [] and loop momenta. The integral has the form 

/^-pW («) 

As emphasized in || there are some very striking features of this result. The first is 
that the pole at P = arises from the high momentum region of integration. Although we 

the rest of the paper, we denote by small p the 4-momentum vector with pi 2 = -Pi. 2 in the 
non-commutativity plane. 
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evaluated it for the massless theory, the pole itself is independent of mass. Furthermore 
this contribution to the self energy has a huge effect on the propagation of long wavelength 
particles. The on-shell condition or dispersion relation becomes 



where c is proportional to the coupling constant. Thus, as discussed in ||, the behavior of 
the non-commutative theory below the non-commutativity scale seems to be nothing like 
the commutative theory. In this case the low momentum end of the spectrum is completely 
removed from the low energy theory. 

Commutative gauge theories are better behaved in the UV than commutative scalar 
theories. The worst divergences in pure Yang mills theory or Yang Mills theory with 
fermions are logarithmic. This naively suggests that in their non-commutative versions 
the worst anomalous effects will be logarithmic in P. As an example consider the vacuum 
polarization correction to the gauge boson propagator. The divergences have the form 



where p 2 is the squared four-momentum of the gauge boson. Note in particular that the 
mass correction vanishes since II vanishes at p 2 = 0. This situation suggests that in the 
noncommutative theory the worst anomalous effect in the propagator has the form 



If this were so, the dispersion relation of a low energy gauge boson would be unaffected 
by the non-commutativity. As we will see in the next section this is generally incorrect. 



In this section we analyze U(l) Yang-Mills theory on a non-commutative space. The 
classical action is given by 





II ~ g 2 p 2 log k 




n ~ g 2 p 2 log P 2 




2 U(l) Non— Commutative Yang— Mills 





with the field strength F given by 



F, 



d^Ay -d u A^- ig[Ay., A v ] 




and 



[^4/^5 A u ] — A^ * A v — A v * A^. 



(2.3) 



3 



The * -product between two functions 4>i(x) and 02 (%) is denned by 



fa * 4> 2 (x) = e^ 9 * a "My)Mz)\y= z = x . (2.4) 

The theory is invariant under non-commutative gauge transformations 

6 x A ll = d ll \-gi(A lt *\-\*A ll ). (2.5) 

We may add matter fields in the theory as well. The scalar and fermionic parts of the 
action that involve interactions of the matter fields with the gauge field are given by 

Smaller = J * -fD^ + \{D^f. (2.6) 

The covariant derivative acts on the fields as follows 

D il X = d ll X-ig[A^X\. (2.7) 

The commutator is defined through the ^-product as before. The matter fields are covari- 
ant under the family of gauge transformations given in Eq(2.5). 

The Feynman rules for the theory have been worked out in references || 1§ . The vertices 
look similar to those of a commutative non-abelian gauge theory with all matter fields in 
the adjoint of the gauge group. The structure constants are replaced by sines of external 
momenta as shown in detail in the Appendix. The Feynman rules for the ghosts are also 
included. The vertices vanish when 9 is taken to be zero as expected. Our computations 
are done in the Feynman gauge. In this gauge the gauge field propagator is given by 

-i 9 ^-. (2.8) 

We refer the reader in [0 for a derivation of the perturbative Feynman rules and gauge 
fixing of the theory. 



3 Photon self— energy correction 

In non-commutative gauge theory the most serious anomalous effects, namely those which 
exhibit inverse powers of P, are associated with the 2 and 3 point functions. We will begin 
with the computation of the 2-point photon self energy diagrams. We will consider the 
contributions from loops involving fermions, scalars and gauge bosons (including ghosts). 
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Figure 2: Photon self energy: gauge loop corrections. 
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Figure 3: Photon self energy: matter loop corrections. 

It should be noted that the individual contributions are gauge invariant. Similar calcula- 
tions have been done by Hayakawa || . The relevant diagrams are shown in Figures 2 and 
3. We will illustrate the procedure with the fermion loop (line (1) in Fig. 3). 
Using the Feynman rules from the Appendix we find 

where Nf is the number of Majorana fermions, each of which counts as two fermion species. 

We are interested in the contribution to the integral coming from very high loop mo- 
mentum. We therefore drop the sub-leading dependence in the integrand and replace 
eq.(3.1) by 

m T (p) = - ig *N f J (^I^P 1 sin* (5*) ■ (3.2) 

Using 

sin 2 (~#) =i[l-cos(pO], (3-3) 

we can isolate the planar and non-planar contributions. The non-planar contribution is 
obtained by dropping the first term and keeping only the cos term. Working out the trace, 
we get 

inj"(p) = tfN, J (2 ^ )4 (2/ / l) (3.4) 
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Since / d A l/ (27r) 4 ^ e^ 1 = ia (log A — log \p\) 0, where A denotes a short-momentum cut-off 
and a is a positive real constant, we can rewrite (|3.4j) as 



iHfip) = 4ig 2 N f a (2& l d u - g» u d 2 ) log \ft\. (3.5) 

The integral Q3.4|) is finite, and the cut-off dependence vanishes after differentiating. We 
note that the term g^ /ft 2 cancels in this expression leaving us with 

mf (p) = -%ig 2 (2N f )a (3.6) 



The answer is somewhat surprising. If the factor sin 2 \\plj were not present in eq(3.1) 
the expression would be the conventional self energy diagram of the commutative theory. 
Gauge invariance would be invoked to say that any quadratic divergence is absent. Alter- 
natively the diagram can be Pauli Villars regulated eliminating the quadratic divergence. 
The integral with the trigonometric factor is finite and well defined. However it quadrati- 
cally diverges as ft — > 0. Thus we see that there is an anomalous effect of order 9~ 2 arising 
out of a diagram which in the commutative theory is quadratically divergent by power 
counting but for which the divergence vanishes as a consequence of symmetry. As noted in 
H in the context of scalar theories this type of behavior is proportional to inverse powers 
of 9. Evidently the limit in which 9 — > does not smoothly tend to the commutative 
theory. 

The physical interpretation of terms like eq(3.6) is very interesting. For small non- 
commutative momentum, the one-loop inverse propagator is given by 



r 



I 2 2 o2\ 2 P/iPf 



(3.7) 



where P represents the projection of the spatial momentum on the (1,2) plane. From 
this matrix we can read the dispersion relation for the two physical, transversely polarized 
photons. Suppose P is along the 2-direction so that P is in the 1-direction. Then the 
photon polarized in the direction perpendicular to P satisfies the same dispersion relation 
as a photon would in the commutative theory 

Vl = Pi + P 2 - (3-8) 

However, the photon polarized along the 1-direction, parallel to P, satisfies a different 
dispersion relation given by 

pl=pl + P 2 + cg 2 -^. (3.9) 



2 In Euclidean space the integral is of course real, / d A l/ {2tt) a ^ e 1 ? 1 = a (log A — log \p\), a > 0. 
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This splitting of the polarization states of the gauge boson is perfectly consistent with 
gauge invariance. Indeed the vacuum polarization tensor in eq.(3.6) is purely transverse 
which follows from the identity pp = 0. This effect would not be possible without the 
breaking of Lorentz invariance caused by 9. 

We remark at this point that since the contributions from the scalar loops and from 
the gauge sector are gauge-invariant by themselves, they individually give combinations 
of the form (|3.4j ). The computations involving gauge bosons (Fig. 2) and scalars (line (2) 
on Fig. 3) in the loop are very similar, and without giving explicit details, we add all three 
sectors together obtaining 



This shows that in the supersymmetric theories with an equal number of bosons and 
fermions, quadratic divergences in the photon self-energy do not appear at one loop. 

As we have mentioned before, the coefficient c = 8(N S + 2 — 2Nf) of the quadratic 
anomalous term is proportional to the number of bosons minus the number of fermions in 
the theory. In particular it is zero in any supersymmetric model. This had to be the case 
since the splitting of the two photon states would be inconsistent with supersymmetry. In 
particular, the fermions in the theory remain massless with dispersion relation 



Note also that in the M = 4 case the SO (6) R-symmetry of the theory insures that the 
scalars do not split. So it would be impossible to achieve bose-fermi degeneracy with 
multiplets of M = 4 supersymmetry had c been non-zero. 

We also note that the coefficient c will vanish in softly broken supersymmetric theories 
as well. In particular, the effect we find arises from high momenta circulating the loop. 
Therefore, it is independent of the mass of the loop particle. As long as the number of 
fermions equals the number of bosons in the theory c will cancel. Simple dimensional 
analysis shows that the only effect of adding a mass term in the Lagrangian is to modify 
the logarithmic singularities as p —>■ 0. 

Finally, we note that c becomes negative if the number of fermions is bigger than the 
number of bosons. This is also the case in the 4 theory when we include fermions coupled 
to the scalar field through Yukawa couplings. In both cases, the theory becomes unstable 
at low energies. 




(3.10) 



j> = 0. 
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Figure 4: Vertex corrections: gauge sector. 

4 Vertex Corrections 

Next we analyze 1-loop corrections to the three-photon vertex. In Figures 4 and 5, we 
have drawn all 1PI Feynman diagrams that contribute corrections to the vertex up to 
1-loop order in perturbation theory. As before, we shall see that there are unfamiliar long 
distance effects arising from the region of integration of high loop momenta. 

Let us consider in detail the scalar graph involving cubic vertices only (Fig. 5, 1st 
graph on line (1)). Applying the Feynman rules, we find 



,3 



dH (2l+p 1 )^(2l-p 2 )^(2l+p 1 -p 2 )^ . fpxl\ . (p 2 l\ . fp^l+px) 



W 7oZ\7 J57i r^TT—^ sm ~7T sm ~TT~ sm 



(2tt) 4 I 2 (l-P2) 2 (l+Pi) 2 \2 J \2 J V 2 / 

(4.1) 

We must also impose overall momentum conservation so that 

Pi+P2+P3 = 0. (4.2) 



We have amputated the external propagators for compactness. Eq (|4.1|) has to be multiplied 
by the number of scalars in the theory. Ignoring the phases for a moment, we note that 
the integrand is linear in I at high loop momentum 

In the commutative non-abelian theory, however, no linear divergence arises because the 
integral is zero by symmetry. Thus it is at most logarithmic in the cutoff. In the non- 
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+ 2 perms + (1) 
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Figure 5: Vertex corrections: matter sector. 



commutative case the oscillating phases spoil the rotational symmetry. They make the 
integral finite but they induce low momentum poles of the form 



pHl pfJ-2 pus 



(4.4) 



Unlike the previous effect which was of order 9~ 2 , this is a Q~ x effect. 

To compute the precise coefficient in front of such anomalous terms, it is enough to con- 
sider high momentum running in the loop ignoring external momenta in the denominators. 
Then the integral becomes 



— 8ig cos' 



P3Pl 



dH %imm>« . (pii\ . (p 2 i\ . fp 3 f 

-. — —r sin — sin — sin — 



(4.5) 



To carry out this integral, it is useful to express the product of sines as a sum of exponen- 
tials. Using Eq(4.2), we can write the product of sines as 



- [sin (pit) + sin (p 2 l) + sin (p 3 /)] . 



(4.6) 



We are left with a sum of three simpler integrals, and, in addition, we can replace each 
sine by an exponential in the integral. In this form, it is easy to see that the contributions 
arise solely from the six non-planar graphs in the double line notation. 

To these, we must add the contributions from the scalar graphs that involve a quartic 
vertex (Second graph and its two permutations on line (1) in Fig. 5). The graphs produce 
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a different tensor structure but the phase factors are the same. There are three such 
graphs from permuting the external particles among the external lines. Each graph has 
to be multiplied by a symmetry factor of 1/2. Adding all four graphs together yields the 
following integrals 

4# 3 iV s cos(^i) / —-7^ UPW 3 - l 2 {l^g^ + perms)! (e tplL + e ip2t + e ml ). (4.7) 
2 J (27r) 4 l b L J 

The pure gauge sector graphs and the fermionic graphs can be computed in the same 
way. One has to remember to include a combinatorics factor of 2 for the ghosts and fermion 
graphs arising from the two different cyclic orderings of the external particles on the loop. 
The gauge boson quartic graph has a symmetry factor of a 1/2. Including these has the 
effect of changing the coefficient in front of the integrals to 

Ns + 2 - 2N f , (4.8) 

where Nf is the number of Majorana fermions in the theory. We see that in any super- 
symmetric theory the linear poles are absent. 

We will now explicitly find the structure of the linear poles. To compute the integral 

in 

" 7 (2vr) 4 Z 6 ' l4 ' yj 

we again use / <i 4 Z/ (27r) 4 ^ e lj5 ' = za (log A — log \p\). Notice that ( |4.y|) is convergent, and 
the A-dependent part can thus be dropped. We can write 

where J(p) = J <i 4 Z/(27r) 4 ^ e Jj5Z = —iJe{p)- The integral over Euclidean space, Je(p) is 
obtained from J(p) after performing a Wick rotation. It satisfies 

n Jb(p) = 4tP 3 p ^G 6 ) = « lo s IpI- ( 4 - n ) 

p a dp dp 

Solving for J and substituting 

AP) = -iJ E {v) = -^P 2 (4 logp - 3) (4.12) 



in ( 4.10p we find 

J)/^1 rjrdi'l qrJ^Z 

jmww = _ 2a y 1 y , (4.13) 



pi 
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Note that the second possible tensor structure, (]5 Ml g^ 2 ^ 3 +p tM2 g^ 1 ^ +p M3 gt 11 ^ 2 ) j P 2 ; cancelled 
completely, just like the photon mass term in the photon self-energy correction ( |3.6|) . This 
fact will turn out to be essential in showing the gauge invariance of the S-matrix in the 
next chapter. 

We can now summarize the computations in this section by giving the linearly divergent 
terms in the correction to the 3-point photon vertex: 



r™(pi,p 2 ,p 3 ) = Sag 3 (N s + 2- 2N f ) cos -p 3Pl £ Pi + . . . , (4.14) 



where the terms denoted by . . . are at most logarithmic in £>i,P2,.P3- 

The anomalous effects found in this paper and in || are highly nonlocal but in a 
particular way. The matrix elements in eqs(3.10) and (4.14) depend depend only on the 
components of momentum in the x 1 , x 2 plane and are independent of x 3 and x A . Thus while 
nonlocal in the noncommutative directions they are completely local in the commutative 
directions. 

In supersymmetric theories the nonlocal 9~ 2 and 9~ x effects are absent, at least at one 
loop. However there are still logarithmic dependences on P which are proportional to the 
one-loop coefficient of the (3 function?] in the corresponding commutative non-abelian gauge 
theory. Once again the corresponding effects are nonlocal only in the noncommutative 
directions. Theories such as N = 4 super Yang Mills theory with vanishing (3 function 
seem to be free of the non-analytic dependence on 9. 



In this section we show that the anomalous terms in the 2-point and 3-point functions 
we computed are consistent with gauge invariance. To check gauge invariance, we study 
a case involving scattering of two gauge bosons into two fermions. We consider 1-loop 
diagrams only and choose the kinematic variables so that the anomalous terms we found 
dominate. 

All particles must be put on shell using the appropriate, corrected dispersion relations 
up to 1-loop order in perturbation theory. We choose one of the gauge bosons to be 
transversely polarized. We denote its polarization by e and its momentum by p; then 



3 The (3 function is the one controlling the running of the t'Hooft coupling in the large TV limit. 




5 Gauge Invariant S— Matrix? 



/* = 0. 



(5.1) 
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Figure 6: Tree level scattering diagrams. 



+ 




Figure 7: The sum of these diagrams is gauge invariant. 

In order to test the gauge invariance we set the polarization of the other gauge boson, 9, 
equal to its momentum q. The momenta of the fermions are denoted by l 2 ). Gauge in- 
variance requires this scattering amplitude to be zero. The tree level diagrams contributing 
to the process are shown in Fig. 6. 

Now we choose a specific kinematic limit. We choose h + h = —I so that I is small. 
Furthermore, we let I 2 be small so that only the s-channel diagrams are important. Two 
types of one loop diagrams are important (Fig. 7), the diagrams involving corrections to 
the intermediate gauge boson's propagator and the diagrams involving corrections to the 
three gauge boson vertex. Corrections to the two fermion-gauge boson vertex are at most 
logarithmic in I and, therefore, sub-leading. We can think of the intermediate gauge boson 
being coupled to some current j M . The exact form of j M is not important for our purposes. 

First consider 1-loop diagrams involving the gauge boson self-energy. In the limit 
when / is small the contribution to the amplitude becomes 



iM 1 = lQtag 4 sin(|)e^ [g^(q - pf + g">(p - J)" + g^(l - qf\ 



(5.2) 



The first factor is a phase factor coming from the 3-photon vertex. Now, using momentum 
conservation — I = p + q, ep = 0, // = and q 2 = 0, we see that this piece becomes 



Z 2 / 4 



(5.3) 



For small /, sin (ql/2j is just ql/2. We note that the dispersion relation of a longitudinally 
polarized photon remains q 2 = since the self-energy corrections we found are transverse. 
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Next we turn to diagrams involving one loop corrections to the vertex. In the limit 
I — > 0, the contribution to the amplitude becomes 

lM2 _ . (5.4, 

This is because in this limit the vertex is dominated by anomalous terms of the form 
\H v l p '/l 4 . Adding the two contributions together, we see that the amplitude becomes zero 
as required by gauge invariance. 

Next we study the case when p — > 0. Again, we study this particular case because we 
can isolate the possible singular terms. The most dangerous 1-loop diagram contributing 
to the scattering amplitude is the s-channel diagram involving corrections to the 3-gauge 
boson vertex. All other diagrams are sub-leading. The 1-loop contribution to the ampli- 
tude is then given by 

iM = -wMMM. (5 . 5) 

We now distinguish between two cases. First we note that if e is perpendicular to p, the 
amplitude is zero. Recall also from the previous section that this is the photon with un- 
corrected dispersion relation. If on the other hand e is along p, the amplitude is not zero 
but given by 

iM, = (5.6) 



This contribution is cancelled by the tree level graph once we use the correct dispersion 
relation for the gauge boson. The tree level graph is given by 

Qp\ (ej)(2 pg) 
2 J P 



^ 3 = -Wsin(^)^^. (5.7) 



The first factor is a phase factor. We also used the following relations: jl = 0, ep = and 
q 2 = 0. Using momentum conservation, we also find that 

2pq = l 2 - p 2 . (5.8) 

Now, when the polarization of the gauge boson is along the p direction, the dispersion 
relation gets corrected as follows 



1 

p 

Therefore, there is order a g A contribution from the tree-level graph given by 

(pj) 
l 2 p 



p 2 = Sag 2 -. (5.9) 



zM 3 = 8iag\qp)^. (5.10) 
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Adding the two together we see that to order g A the amplitude is zero as it is required by 
gauge invariance. 

The case involving scattering of scalars is more subtle because the two scalar - gauge 
boson vertex also contains linear divergences in 9. We defer this case for an upcoming 
paper 0. 

6 Conclusions 

The most naive expectation about the non-commutative field theories is that they become 
commutative when the non-commutativity parameter 9 % i —>■ limit is taken. In other 
words, star-products are replaced by ordinary products of the fields in the Lagrangian in 
this limit, and one may expect that the theory becomes commutative also at the quantum 
level. As it was shown in || this is generally not true due to the appearance of the 
new divergences at low non-commutative momenta. In the commutative theories, these 
divergences appear at high momenta in the superficially divergent loop integrals, but they 
can be eliminated by an appropriate choice of the regularization scheme. In the non- 
commutative theories, some of these divergences simply do not occur due to oscillating 
phases associated with star-products in the vertices which make the integrals finite. The 
non-commutative momenta thus play the role of the regulator. The dependence on these 
non-commutative momenta, however, does not disappear and manifests itself in the form 
of the new infrared divergences at small values of non-commutative momenta. These 
effects can also be characterized as non-analytic behavior in 9. 

A less naive expectation would be that a non-commutative gauge theory must be free 
of quadratic and linear poles at low non-commutative momenta since the corresponding 
commutative non-abelian gauge theory contains at most logarithmic divergences. In this 
paper we have shown that even this expectation does not hold, and both quadratic and 
linear poles appear in a generic gauge theory. The structure of these new poles is con- 
sistent with gauge invariance but not Lorentz invariance. These effects are local in the 
direction perpendicular to the non-commutativity plane and completely non-local in the 
non-commutative directions. 

In supersymmetric gauge theories these poles cancel between the bosons and the 
fermions at the one loop level but even these theories typically contain logarithmic di- 
vergences at small values of non-commutative momenta. We expect that in AT = 4 SYM 
theory even the logarithmic divergences do not occur, and thus this theory is completely 
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free of anomalous effects in the small non-commutative momentum limit. This is the only 
theory that we know that reduces to its commutative counterpart in the limit y — > 0. 
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-2g sin(|pip 2 )x 

[(Pi - P2^ 3 9^ 2 + (P2 - P 3 ) m ^ 2At3 + (P3 - Pl)^ 3 " 1 ] 



sin(|pip 2 ) sm(|p3j9 4 )(^ 1 ' 13 ^ 2 ' 14 - j" 1,,4 f " 3 )+ 
sin(±p 3 pi) sin(|p 2 P4)(^ 1 ' 14 ^ 2 ' 13 - g^^g^ 4 ) + 



sin(-piP2) 



2tf(Pi +P 2 t )sin(-p 1 p 2 ) 



x sin(ipip 3 ) sin(|p 2 p 4 ) + sin(±pip 4 ) sin(|j5 2 p 3 ) 



2^ sin(-piP2) 
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